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1. Introduction 

There has recently been interest in studying quantum field theories with a nontrivial in- 
frared fixed point (IRFP), both in continuum and on the lattice. Under the renormalization 
group evolution, the coupling of these theories shows asymptotic freedom at small distances, 
analogously to QCD, but flows to a fixed point at large distances where the theory looks 
conformal. Such theories have applications in beyond Standard Model model building. 
These include unparticles, i.e. an infrared conformal sector coupled weakly to the Stan- 
dard Model ||, ||, H|, H|, and extended technicolor scenarios, that explain the masses of the 
Standard Model gauge bosons and fermions via strong coupling gauge theory dynamics 
11 0,1- 

In addition to direct applications to particle phenomenology, the phase diagrams of 
gauge theories, as a function of the number of colours, N, flavours Nf and fermion rep- 
resentations, are interesting from the purely theoretical viewpoint of understanding the 
nonperturbative gauge theory dynamics from first principles. In figure |l| we show a sketch 
of such phase diagram for SU(iV) gauge theory ||. In addition to the fundamental repre- 
sentation, the figure shows the phase structure in the cases of two-index (anti)symmetric 
and adjoint representations. The shaded regions in the figure depict the conformal win- 
dows for each of these fermion representations; below the conformal window the theory is in 
the chiral symmetry breaking and confining phase, while above the conformal window the 
theory is in the non-Abelian QED-like Coulomb phase. The upper boundary in each case 
corresponds to the loss of asymptotic freedom, i.e. to the value of Nf where the one-loop 
coefficient /3q = U/3N — 4/3NfT(R) of the /3-function vanishes. The group theory factor 
T(R) is defined for each representation as 

Tr(T a T b ) = T(R)5 ab . (1.1) 
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Figure 1: The phase diagram of SU(iV) gauge theory as a function of the number of colours, 
flavours and fermion representations (F = Fundamental, 2A = 2-index antisymmetric, 2S = 2- 
index symmetric, Adj = Adjoint). The shaded bands indicate the estimated conformal windows. 



For fundamental (F), two-index symmetric (2S), two-index antisymmetric (2AS) and ad- 
joint (A), the concrete values are, respectively, T(F) = 1/2, T(2S) = (N-2)/2, T(2AS) = 
{N + 2)/2 and T(A) = N. 

For values of N and Nf inside the conformal window, the theory is expected to have 
a non-trivial IRFP. Near the upper boundary one expects the value of the coupling at the 



fixed point, a* , to be small and perturbation theory to be applicable [10]. However, as one 
moves deeper into the conformal window, by e.g. lowering Nf at fixed N, the fixed point 
coupling grows, and nonperturbative methods are required for the analysis. In figure H we 
have used the traditional estimate of evaluating the critical coupling for chiral symmetry 
breaking in the ladder approximation and setting it equal to the fixed point value of the 
two-loop coupling. In other words, at the lower boundary of the conformal window the 
fixed point coupling becomes super critical with respect to the critical coupling for the 
onset of chiral symmetry breaking. 
Explicitly, at two loops we have 



A) = 
Pi = 

which implies that IRFP is at 
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Pi 



(47T). 



(1.2) 



(1.3) 



Above, C*2(i?) is the quadratic Casimir operator in representation R and, for the representa- 
tions we consider, it assumes values C 2 (F) = (N 2 -1)/(2N), C 2 (2(A)S) = (N±1)(N^2)/N 
and C 2 (A) = N. The estimated critical coupling on the other hand is a c = n / (3C 2 (R)) 
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1 11]. Solving a* = a c , one obtains the lower boundary of the conformal window, denoted 
by dashed lines in figure |. The ph ase diagram in the figure |] is therefore a conjecture 
which must be checked by nonperturbative analysis. Since currently the lattice simulations 
provide the only robust nonperturbative method for non-super symmetric four dimensional 
gauge theories, such a check provides an interesting challenge for the lattice community. 

In this work we study SU(2) gauge field theory with Nf = 4, 6 and 10 massless flavours 
of fermions in the fundamental representation on the lattice using 0(a) improved Wilson- 
clover fermions. On the basis of two-loop beta-function six and ten flavour SU(2) theories 
may have a fixed point and be within the conformal window, while the four flavour theory 
is expected to confine. The model with ten fermions is chosen because it is close to the 
upper edge of the conformal window (Nf = 11, where asymptotic freedom is lost) and is 
expected to have a fixed point at a rather small coupling. Thus, in this case we can compare 
the results from lattice simulations with those obtained from perturbative computations. 
The universal 2-loop /3-function has a zero at g 2 ~ 2.90, and at 3 or 4 loops in the MS- 
scheme it is at g 2 ~ 2.47 or 2.52, respectively jL2|. While the MS-scheme results cannot be 
directly compared with the Schrodinger functional results, the convergence indicates that 
the perturbative result should be reasonably accurate. 

In the model with six fermions, however, the fixed point determined by the two loop 
beta function is at much larger coupling, where we cannot expect the perturbation theory 
to hold. Indeed, the location of the fixed point varies significantly at different loop levels: 
using the MS-scheme /3-function at 2, 3 and 4 loops, the location is <?| p « 140, 21 and 
30, respectively. Such a large value of the fixed point coupling may be super critical with 
respect to the critical coupling for the onset of spontaneous chiral symmetry breaking; this 
is explicitly illustrated in figure [l], where the two colour and six flavour theory is already 
outside the conformal window. However, due to the nonperturbative nature of the lower 
boundary of the conformal window, one cannot exclude the possibility that the six flavour 
theory is very close or even within the conformal window. In addition, if we treat Nf 
as a continuously variable parameter, the fixed point vanishes from the 4-loop /3-function 
when Nf is lowered only slightly from 6 to 5.945. In the proximity of this value the 4- 
loop /3-function obtains a shape typical of walking coupling, i.e. the /3-function has a local 
maximum with a small negative value. Clearly, these results imply that the six flavour case 
can be resolved only with non-perturbative calculations. 

On the other hand, the four flavour case is expected to be QCD-like, with confinement 
and chiral symmetry breaking. It is included here for comparison with the larger Nf cases. 

On the lattice, SU(2) gauge field theory with Nf > 2 fundamental representation 
fermions has been recently studied in Refs. [13, 14, 15 \} The study by Bursa et al. jD| of 



SU(2) with six fundamental flavours suggested a possibility of an IRFP at much smaller 
coupling than expected from perturbation theory. However, this work used unimproved 
Wilson fermions which can be expected to be subject to large discretization errors. Indeed, 



1 In related work, the existence of the infrared fixed point in SU(2) gauge theory with two adjoint 
representation fermions (which is of interest for technicolor model building) has been studied in [fujl [l?], 



§. §. 1 

flavours in 



24, |25|, |2q]. and SU(3) gauge with various fermion representations and numbers of 

|], § || H 0, || H H EH H H- 



-3- 



as described below, the leading order perturbative analysis reveals that the running cou- 
pling measurement (step scaling using Schrodinger functional) using unimproved Wilson 
action can have finite cutoff effects of order 30-60%, whereas the improved action reduces 
these to a few percent level using lattices presently within computational reach. In our 
simulations using 0{a) improved Wilson-clover fermions we indeed observe large deviations 
from the unimproved results. 

We measure the running coupling using the Schrodinger functional method in theories 
with four, six and ten flavours, and the mass anomalous dimension in the phenomenolog- 
ically most interesting six flavour case. Unfortunately, in the six flavour case we are not 
able to fully resolve whether the fixed point exists, but the possible locations of the fixed 
point moves to significantly larger coupling (g 2 J> 10) than the unimproved lattice action 
results indicate. 

This paper is structured as follows: In section ^ we introduce the model and describe 
the Schrodinger functional method, and present the perturbative step scaling function 
results. The lattice simulations and results are presented in section || and in section ||| we 
conclude. 

2. The model and the Schrodinger functional method 

The theory is defined by the action 

S = S G + S F , (2.1) 

where Sq is the standard Wilson single plaquette gauge action for the SU(2) Yang Mills 
theory 

S G = Pl f 1 " \ Tl P»(x)U v {x + ap)Ul(x + av)Ut(x)]\ , (2.2) 
with fa = The part Sf is the clover improved Wilson fermion action 



— — % 

ip a (x)(iD + m )ip a (x) + ac sw tp a (x)-a IJ , u F^ u (x)ip a (x) 



a=l x 

where D is the standard Wilson-Dirac operator 



(2.3) 



^ = ^MV; + V M )-aV*V M ], (2.4) 

with the gauge covariant forward and backward lattice derivatives 

V^(x) = l/a[U^{x)i){x+afL)-il>{x)l V*^(x) = l/a[^{x)-Ul{x-afi)^{x-afi)}. (2.5) 

The clover improvement term contains the usual symmetrized field strength tensor. We 
set the improvement coefficient c sw to the iVy-independent perturbative value p3| 

tW = l + O.1551(l)gg+O(0$). (2.6) 
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While c sw can be determined nonperturbatively 1 44] , in our tests we observed that at strong 
lattice coupling the nonperturbatively determined c sw is close to the perturbative one at 
Nj = 6 and 10. The situation is very different at Nj = 2 [45, where c sw diverges as g$ 



increases. Thus, while the perturbative result ( |2.6| ) does not give full cancellation of the 
0(a) discretization effects, it is a concrete recipe which we expect to eliminate most of the 
O(a) effects. We also include the perturbative improvement at the Schrodinger functional 
boundaries as described in [47, 46 1. 



We measure the running coupling using the Schrodinger functional method [48, 49, 50 



51] . The coupling is defined as a response to the change of the background field and the 
scale is set by the finite size of the lattice. We consider a lattice of volume V = L 4 = (Na) 4 . 



The spatial links at the t = and t = L boundaries are fixed to [49 

U^x^t = 0) = e~ i ^ a / L (2.7) 

U^X, t = L)= e -^~r,)a 3 a/L 

with CJ3 the third Pauli matrix. The spatial boundary conditions are periodic for the gauge 
field. The fermion fields are set to vanish at the t = and t = L boundaries and have 
twisted periodic boundary conditions to spatial directions: ip(x + Li) = exp(i7r/5)^(a;). 

At the classical level the boundary conditions generate a constant chromoelectric field 
and the derivative of the action with respect to rj can be easily calculated: 

drj g£ 



where k is a function of N = L/a and rj [49|. At the full quantum level the coupling is 
defined by 

W = T < 2 ' 10) 

To quantify the running of the coupling we use the step scaling function H(u,s,L/a) 
introduced in pq| . It characterizes the change of the measured coupling when the linear 
size of the system is changed from L to sL while keeping the bare coupling (and hence 
lattice spacing) constant: 

E(u,s,L/a) = 9 2 {go,sL/a)\ g2[g ^ L/a)=u (2.11) 
a(u,s)= lim E(u,a,L/a) (2.12) 

a/L->0 

In this work we choose s = 2. To obtain the continuum limit cr(u, s) we calculate 
T,(u, s, L/a) at L/a = 6 and 8. Since we expect the discretization errors to be (mostly) 
removed to the first order in the improved action, we use quadratic extrapolation to find 
the limit a — > 0. 

The step scaling function is related to the /3-function by 

- 21nW = /„ Trnvry (2 ' 13) 
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Close to the fixed point, where the running is slow and \f3\ small, we can approximate the 
/3-function by 



2m(s) 



9 2 



(2.14) 



The estimating function (3*(g) is exact at a fixed point but deviates from the actual /3- 
function as \g — g*\ becomes large. 

We also measure the mass anomalous dimension 7 = d\nm q / d\n /U of the theory with 
6 fermion flavours using the pseudoscalar density renormalization constant. The mass 
anomalous dimension is of interest for this case, because the possible infrared fixed point 
is at non-perturbative value of the coupling. It is defined on the lattice as [52| 



ZpiL) - ML/2)' 
where /1 and fp are correlation functions of the pseudoscalar density 

-1 



h 
fp(xo) 



12L 6 
-1 
12L 6 



d\d 3 vd 3 yd 3 z(C(uh 5 \ a C(vK(yh 5 \ a ((z)), 
d 3 yd 3 z(^(x )75AV(^o)C(y)75A a C(z)). 



(2.15) 

(2.16) 
(2.17) 



For these measurements the boundary matrices at t = and t = L are set to unity; thus, 
separate simulations are needed. The mass step scaling function is then defined as in p3|: 



Ep(n, s, L/a) 



Z P (g ,sL/a) 



Zp(go,L/a) 



9 2 {90,L/a)=u 



ap(u,s)= lim T,p(u, s, L/a), 

a/L->0 



(2.18) 
(2.19) 



and we choose again s = 2. We find the continuum step scaling function ap by measuring 
Sp at L/a = 6, 8 and 10 and doing a quadratic extrapolation. 

The mass step scaling function is related to the anomalous dimension by (see [152]) 



a P {u,s) 



a(u, s) 



do/(26 ) 



exp 



- dx 


( j(x) 


_ M 


Jy/u 




b x J 



(2.20) 



where 60 = A)/(167r 2 ) in terms of the one- loop coefficient [3q = 22/3 — 2iVy/3 of the beta 
function and do = 3C2(F)g 2 /(8ir 2 ) = 9/(32-7r 2 ) is the corresponding one-loop coefficient 
for the anomalous dimension, 7 = —dog 2 . Close to the fixed point ( |2.20| ) simplifies con- 
siderably. Denoting the function estimating the anomalous dimension j(u) by J*(u), we 
have 



log crp{g 2 ,s) = -j*(g 2 ) 
^i*{g 2 )-- 



1 



-7*(ff 2 )logs, 

M 

log a P (g 2 ,s) 
log s 



(2.21) 
(2.22) 
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Table 1: Parameter k used in the simulations at each /3l = 4/<?o an d the number of measurements 
performed on the largest lattice. 



The estimator 7*(g 2 ) is exact only at a fixed point where fi(g 2 ) vanishes and deviates from 
the actual anomalous dimension when (3 is large. 

The Schrodinger functional boundary conditions enable one to run at vanishing quark 
mass. The Wilson fermion action breaks the chiral symmetry explicitly and allows additive 
renormalization of the quark mass. We therefore find the quark mass from the PCAC 
relation 

M( , _ l (d* +d )f A (x ) a d* d f P (x ) 
aM(X0) "4 /p(* ) +CA 2 /p(x ) • ^ 

We have used here the improved axial current 

fA = fA + ac A ±(d; + d IM )f P , (2.24) 
fA(x ) = Y^j d 3 yd 3 z(^(x ) lol5 X a ^(x )C(yhol 5 X a C(z)). (2.25) 

In this work use the perturbative expansion for the improvement coefficient c A [43]: 



CA = -0.00567(l)C 2 (F)gg + 0(g%), C 2 (F) = 3/4. (2.26) 

The Schrodinger functional boundary conditions remove the zero modes that would 
normally make it impossible to run simulations at zero mass. We define n c as the value of 
the parameter k where the mass aM(L/2) vanishes. To find k c we measure the mass at 
3 to 7 values of k on lattices of size L/a = 16 and interpolate to find the point where the 
mass becomes zero. We then use the same value of k c on all lattice sizes. The values of k c 
used in the simulations are given in table EL In practice we achieve |aM| < 0.01. 
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Figure 2: Contribution of a massless Wilson quark to the step scaling function normalized to its 
continuum value at one loop order in perturbation theory. The top three curves show the result 
for gauge groups SU(4), SU(3) and SU(2) (from top to bottom) for unimproved Wilson fermions, 
while the lower three curves show the result after 0(a) improvement has been taken into account. 



3. Improvement in perturbative analysis 

The degree of improvement obtained with the clover term and boundary counter terms 
can be quantified with one loop order perturbative analysis of the step scaling function 
(2.11). Because the gauge action is identical with improved and unimproved fermions, to 



this order it is sufficient to consider only the fermion contribution to the step scaling. To 



one loop order the step scaling function (|2.1l[) can be expanded as 



E(w,s,L/a) = g 2 {go,sL/a)\ g 2 {gihL/a)=u 

= n + (S li0 + S 1 , 1 iV / )n 2 . (3.1) 

The fermion contribution is denoted by En. To evaluate these perturbative contributions 
we use the methods in [p4| , |55| , and choose s = 2. The continuum limit of E^i is given by 
the fermionic contribution to the one loop coefficient bo = (3q/(1Qtt 2 ) of the beta function, 
i.e. 

lim En = 2& ,iln2, (3.2) 

Lja— >0 

where 6 ,i = l/(24vr 2 ). 

In figure ^ we show 5 = Si i x/(26o,i In 2) both for unimproved Wilson fermions and 
with 0{a) improvement. One immediately observes that without improvement, E^i de- 
pends strongly on L/a and approaches the continuum limit only for presently impractically 
large lattices, while with improvement the large lattice artefacts are absent. While this 
level of improvement is probably not preserved at higher orders, this nevertheless strongly 
motivates the use improved actions in the lattice studies of these theories with Wilson 
fermions. 
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Table 2: The measured values of g 2 at each Pl — ^/g 2 an d L/a with 4 flavours of fcrmions. 
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Table 3: The measured values of g 2 at each /?£ and L/a with 6 flavours of fermions. 



4. Simulations 

We use hybrid Monte Carlo (HMC) simulation algorithm with 2nd order Omelyan inte- 



grator [56, 57] and chronological initial values for the fermion matrix inversions p8| . The 
trajectory length is 1, and the step length is tuned to have acceptance rate larger than 
80%. 

The measured values of the running coupling squared in theories with 4, 6 and 10 
fundamental representation fermions are given, respectively, in tables ||, |3| and ^ and shown 
in figures |3| and The measurements are taken after every trajectory, and the number 
of trajectories for each point is up to 230,000 (for Nf = 6, volume 16 4 and largest 
values). For 6 and 10 flavours the strongest lattice couplings (smallest /3^) are very close to 
the strongest practical values for our action; at smaller /3l the simulations become either 
too slow or unstable, possibly signalling a proximity of a bulk phase transition. These 
transitions are a lattice artefact and limit the range of allowed lattice couplings. 
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Table 4: The measured values of g 2 at each (3l and L/a with 10 flavours of fermions. 




Figure 3: The measured values of g 2 (go, L/a) against a/L with 4 and 6 flavours of fermions. The 
black dashed line gives an example of the running in 2-loop perturbation theory at modest coupling, 
normalized so that it matches the measurement at L/a = 6. 



One can immediately recognize the main features from figures [| and |j: at Nf = 4, 
the coupling becomes stronger as lattice size increases, and the running becomes faster at 
stronger coupling. This agrees with the expected QCD-like behaviour. In contrast, with 10 
flavours we observe basically no running at all within our statistical errors, except with the 
strongest lattice coupling = 1) used, where g 2 becomes smaller as volume increases. 
This is the expected behaviour if we are at the strong coupling side of an infrared fixed 
point, where the /3-function is positive. However, at /3l = 1 the data differ qualitatively 
from the larger /3l measurements, possibly a signal of contamination from finite lattice 
spacing effects at large bare coupling. 

In the theory with 6 fermions the running remains small and slightly positive in the 
studied range of couplings. To illustrate the running we show the scaled step scaling 
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Figure 4: The measured values of g 2 (<?o, L/a) against a/L with 10 flavours of fermions. The black 
dashed line gives an example of the running in 2-loop perturbation theory. 
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Polynomial 
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6 


Polynomial 
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30.20 


16.11 
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Rational 
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7 


3.545 


2 
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Polynomial 


13.744 
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13.802 


5.732 
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Rational 


9.896 


6.271 


12.704 


6.301 
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Table 5: \ values and free degrees of freedom in the polynomial fits 4.1 and the rational fits 4.2 



function E(g 2 ,2,L/a)/g 2 = g 2 (g 2 ,2L/a)/g 2 (g 2 , L/a) in figure | at L/a = 6 and 8. The 
running is compatible with the perturbation theory at small g 2 , but deviates from it at large 
coupling. The scaled step scaling measurements decrease at g 2 > 4 and may reach unity at 
g 2 > 10, indicating an ultraviolet fixed point. However, there appears to be a systematic 
difference between the L/a = 6 and 8 points at large g 2 . Thus, proper continuum limit is 
necessary. 

We also note that at L/a = 6 the very largest coupling (3l = 4/g 2 , = 1.39 point 
deviates substantially from other points. This is caused by the finite size effects at the 
smallest volume, clearly visible in figure ||. 

The continuum limit of the step scaling function T,(g 2 ,2, L/a) has to be evaluated at 
constant Schrodinger functional coupling g 2 . However, the measurements were performed 
at selected fixed values of g 2 . = 4//3l, which do not correspond to the same <7 2 -values at 
L/a = 6 and 8. Therefore, it is necessary to shift the measurements so that (/ 2 -vaiues 
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Figure 5: The Nf = 6 scaled lattice step scaling function 'S(g 2 ,2,L/a)/g 2 = 
g 2 (gQ, 2L/a)/g 2 (g 2 ,, L/a) for Nf = 6 calculated directly from the data in table ||. The black dashed 
line is the continuum 2-loop perturbative result for a(g 2 , 2). 
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Figure 6: The interpolating functions for six fermions. The plot on the left shows the polynomial 
interpolation and the plot of the right shows the rational interpolation. For clarity, the graphs 
are displaced up by 3,2,1 and and left by 0.3,0.2,0.1 and for lattice sizes L/a = 6,8,12 and 16 
respectively. 



match. This could be done by performing new simulations so that g 2 at L/a = 6 and 8 
match, or by reweighting in and k. However, a much more economical and convenient 
way to achieve this is to interpolate the actual measurements of g 2 (gQ 1 L/a) at each lattice 
size L/a by fitting to a function of <?q. This results in a "measurement" of g 2 (<7g, L/a) over 
a continuous range of g\ . 



-12- 



It is necessary that the measurements cover the interpolated range densely enough; 
otherwise the form of the fitted function ansatz may influence the final results. Obvi- 
ously, the fit must also be statistically good in order to describe the true behaviour of the 
data. The interpolating function also somewhat averages out fluctuations in individual 
data points. 

We have done the interpolation using polynomial and rational interpolating functions. 
Polynomial function is a power series in the bare coupling, and it is commonly used in the 



literature |2|, |§: 



1 1 



g 2 (g 2 l ,L/a) gl 



(4.1) 



i=l 

In the case of 4 and 10 flavours we truncate the series at n = 4 and in the case of 6 flavours, 
where more /^-values are available, at n = 5. These values were chosen to optimise the 
confidence levels of the fits, while keeping the number of the fit parameters tolerably small. 
We also kept the same number of terms for each L/a. The x 2 -values of the fits are shown 
in table |5|, and, for Nf = 6, the fit is shown in figure ||, left. 

It is clear that the polynomial fits do not fit the data very well, especially at the 
strongest bare couplings where the measured g 2 increases very rapidly. This leads us to 
try a rational function interpolation: 



1 _ 1 

g 2 (gl,L/a) gl 



For 4 and 6 fermions the number of terms were chosen to be n = m = 2 and for 10 fermions 
n = 1 and m = 2. 2 Again, these parameters were chosen to optimize the confidence level 
without producing singularities within the range of gl of the measurements. The stability 
of the fits was checked by varying n or m. The rational interpolation captures the rapid 
increase at strong bare couplings better, and the confidence levels of the fits are better in 
all cases, as shown in table || and figure ||. Thus, we perform the subsequent analysis using 
the rational interpolation functions. 

The interpolating functions are used to calculate s,L/a) with L/a = 6, 8 in a 
continuous range of values of u = g 2 . This enables us to perform a continuum extrapolation. 
Since we expect most of the order a errors to have been removed, we fit the data at L/a = 6 
and L/a = 8 with a function of the form 

£(u, 2, L/a) = a(u, 2) + c(u) (L/ a y 2 (4.3) 

independently at each value of u. This extrapolation is shown for Nf = 6 for three se- 
lected values of u in figure |7|. To estimate the systematic errors from the extrapolation 
we compare with the result obtained by simply taking the largest volume step scaling, 
a(u,2) = S(n,2,8). 

Because we only have step scaling data at L/a = 6 and 8, we are naturally not able to 
verify that the first order in L/a -term is indeed small or that the higher order contributions 



2 For Nf — 10 choosing n — 2 produces a singularity within the fitting range, which is not acceptable. 
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Figure 7: The step scaling function extrapolated to the continuum limit with Nf = 6 and for 
three chosen values of u = g 2 . 



can be neglected. However, we note that using a first order only extrapolation would give 
results roughly comparable with the extrapolation (|4.3|), but with larger statistical errors 
and further away from the largest volume result. The L/a = 10 results with Nf = 6 are 
not used for the measurement of the step scaling for the running coupling, but are needed 
for the measurement of the anomalous dimension. 

The final results of the step scaling functions are shown in figures g | and for 
4, 6 and 10 fermion flavours respectively. The errors shown include only the statistical 
errors from the measurements, fits and extrapolation. A rough measure of the systematic 
errors can be obtained from the comparison of the infinite volume extrapolated result and 
the result using only the largest volume step scaling without extrapolation. The error 
propagation is calculated using jackknife blocking throughout the whole analysis. 

We have also calculated the estimator for the ft- function, ft*, defined in Eq. ( 2.14| ). 
These are shown in figure 11, together with the 2-loop ft- function and ft* obtained from 
the 2-loop perturbative step scaling function. Thus, the difference between the perturbative 
curves gives a measure of the error made in the approximation in Eq. (2.14). 

What can we conclude from the results? From figure [8| we can observe that, in the 
range of couplings studied, the Nf = 4 case behaves as expected: /3-function is negative 
and becomes smaller as the coupling increases, and agrees overall with the 2-loop perturba- 
tive /3-function. Naturally deviations from the perturbation theory are expected at larger 
couplings. However, for our purposes it is sufficient to verify the QCD-like behaviour and 
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Figure 8: (colour online) The scaled step scaling function a(g 2 ,2)/g 2 with 4 fermions. The thick 
red line corresponds to the continuum extrapolation, Eq. (|4.3|), and the hashed band to the statistical 
errors of the extrapolation. The thick dashed line with the shaded error band is the largest volume 
step scaling function without extrapolation. The thin dashed line is the 2-loop perturbative value 
ofa(g 2 ,2)/ 3 2 . 
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Figure 9: As in figure @ but with 6 fermion flavours. 



we did not attempt to reach smaller couplings. 

For Nf = 6, the results in figure ^ show that the largest volume step scaling deviates 
from the perturbative result already at g 2 ~ 4, after which it decreases. However, the 
continuum extrapolation remains close to the perturbative curve up to g 2 ~ 7-8. After 
this the curve turns sharply downwards; this behaviour is caused by the anomalous L/a = 6, 
Pl = 1.39 point visible in figures || and [|. The largest volume step scaling (shaded band) 
can be compared to the corresponding L/a = 8 uninterpolated step scaling in figure ||. 
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Figure 10: As in figure || but with 10 fermion flavours. 



The large difference between the largest volume and extrapolated step scaling functions 
is naturally due to the fact that the step scaling T,(g 2 , L/a) at L/a = 6 and L/a = 8 differ 
substantially, and the lever arm to the continuum {a/L — > 0) is long. The variation 
between these two curves gives an estimate of the systematic errors in the extrapolation. 
The largest values of the Schrodinger functional coupling we achieved was around g 2 ~ 
12 (a ~ 0.96). Unfortunately, neither the coupling is strong enough nor the errors are 
sufficiently small in order to unambiguously distinguish between an infrared fixed point at 
g 2 > 11 and "walking" behaviour, where the step scaling function starts to increase again 
at stronger coupling. With our choice of the action we cannot reach stronger couplings 
because simulations become rapidly impractical at smaller fii than 1.39, our strongest 
lattice coupling at Nj = 6. Obviously, calculating the step scaling at volumes L/a = 10 
and 20 (or larger) would be needed in order to stabilize the continuum extrapolation. 

Finally, for Nf = 10 the measured /3-function is essentially compatible with zero or with 
the perturbative /3-function at g 2 < 2.5. In this case the 2- loop perturbative /3-function 
is expected to be relatively accurate (higher order perturbative corrections in MS-scheme 
are small (l^]). The value of the /3-function is very small, and clearly, the accuracy of 
our results falls far short from being able to resolve the behaviour of the /3-function in 
this range. Nevertheless, theoretically we know that the /3-function must be negative at 
small coupling. At g 2 j> 2.5 the measured /3-function deviates significantly upwards from 
the perturbative one. Combined with the knowledge that the /3-function is negative at 
small couplings, this indicates the presence of an infrared fixed point. The onset for this 
behaviour is very close to the fixed point in the 2-loop /3-function. This deviation is in 
practice caused by the strongest coupling (3^ = 1 data from the lattice: from figure^] we can 
see that only at (3l = 1 the measured coupling g 2 > 2.5. At these couplings we can expect 
strong lattice artefacts; thus, we believe that the deviation from the perturbative value 
is caused by unaccounted for systematic errors in the continuum extrapolation. We still 
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Figure 11: The estimators for the j3- function in equation Eq. ( |2.14 ) 
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L/a = 10 
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L/a = 16 


L/a = 20 


2.4 


0.9666(8) 


0.9353(9) 


0.9171(19) 


0.9014(19) 


0.8870(15) 


0.865(4) 


2 


0.8953(10) 


0.857(3) 


0.838(2) 


0.823(3) 


0.793(4) 


0.766(6) 


1.5 


0.702(4) 


0.669(3) 


0.646(4) 


0.618(6) 


0.586(5) 


0.573(6) 


1.44 


0.636(3) 


0.610(3) 


0.588(4) 


0.572(4) 


0.548(4) 


0.517(6) 


1.4 


0.543(5) 


0.547(5) 


0.539(5) 


0.534(5) 


0.508(6) 


0.480(8) 


1.39 


0.508(5) 


0.515(7) 


0.520(6) 


0.517(6) 


0.488(8) 


0.476(9) 



Table 6: The measured values of Zp at each /3l and L with 6 flavours of fermions. 



emphasize that in absolute numbers the deviation is still rather small; in order to resolve 
very slow running we would need extremely accurate measurements with correspondingly 
very small systematic errors. This does not appear to be doable at Nf = 10 with the 
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L/a = 6 


L/a = 8 


L/a = 10 


L/a = 12 


L/a = 16 


L/a = 20 


25.15 


6.48 


0.28 


0.84 


3.00 


1.50 



Table 7: x values (2 degrees of freedom) for the Zp fits for Nf — 6 using Eq. (4.4) 



methods used here. 

We measure the anomalous dimension for the interesting case of SU(2) with 6 flavours 
of fermions using the mass step scaling method described in section 0. This method is much 
less noisy than the measurement of the coupling, and we can now reliably use lattices of 
size 20 4 . The measurements of Z p are listed in table |6[ 

We start the analysis by finding an interpolating function to the measured values of 
Zp by fitting with a power series 

n 

Z P (f3 L ,L/a) = l + Y J Ci9 2 \ (4.4) 
i=i 

where we have truncated the series at n = 4. The x 2 values for the fits are given in table 
|7[ The fit at L/a = 6 has bad x 2 / value; however, its effect to the final extrapolation is 
small. We check the systematic errors of the interpolation by also truncating at n = 3; 
the results remain essentially the same with somewhat increased statistical errors. From 
the interpolated Z p ((3l, L/a) we obtain the mass step scaling function £p(u, s, L/a) at 
L/a = 6,8 and 10 using Eq. ( |2.18| ), using u = g 2 from rational fit in Eq. ( ^4.2| ). The 
continuum extrapolation is then done by fitting to the extrapolating function 

T,p(u,2, L/a) = crp(u, 2)+c(u)(L/ay 2 . (4.5) 

The fit is shown in figure We check the systematic effects in the extrapolation by 
leaving out the smallest volume step scaling at L/a = 6. The results remain compatible 
with each other, but with naturally much smaller final statistical errors when L/a = 6 is 
included. Step scaling function is converted to the estimate of the anomalous dimension 
using Eq. ( |2~2~2l) , and the results are shown in figure KL3L 



In this case the non-trivial continuum extrapolation is essential in order to find the 
monotonic growth of 7*(<? 2 ). If we would take e.g. L/a = 8 step scaling alone to estimate 
the continuum behaviour, 7* would decrease at g 2 > 6. This emphasizes the importance 
of the controlled continuum limit. 

The mass anomalous dimension we measure is somewhat smaller than the perturbative 
one at strong coupling. It remains small at all measured values of the coupling, and if there 
is a fixed point at g 2 > 12, the anomalous dimension is > 0.25. 

The behaviour of the anomalous dimension is in contrast with the results obtained 



with an unimproved Wilson action [13], where 7* significantly larger than the perturbative 
one was seen already at g 2 ~ 5, albeit with large errors. Both of the measurements were 
done using the Schrodinger functional scheme, and thus should yield identical results in the 
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Figure 12: The mass step scaling function extrapolated to the continuum limit, using Nt = 6 
data and shown for three chosen values of g 2 . 
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Figure 13: The estimate of the mass anomalous dimension, j*(g 2 ), with Nt = 6. Shown are the 
continuum extrapolation using the mass step scaling S(L/a, 2,g 2 ) with all three volumes L/a = 
6, 8, 10 (hashed band) and only with L/a — 8, 10 (shaded band). 



continuum limit. This again points to the importance of controlling the cutoff effects and 
the continuum extrapolation. On the other hand, the behaviour we measure is qualitatively 
similar to the one observed in SU(2) gauge theory with adjoint fermions using HYP-smeared 
improved fermions [26|. 

For the Nj = 10 case we did not perform a measurement of the mass anomalous 
dimension. Since our results for the coupling constant in the Nf = 10 case are compatible 
with perturbative results, we expect this to hold also for the mass anomalous dimension. 
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In other words, we expect that 7* = 0.08 corresponding with a perturbative fixed point at 
g 2 ~ 2.90. 

5. Conclusions 

In this paper we have presented the results of a lattice study of SU(2) gauge theory with 
Nf = 4, 6 or 10 flavours of Wilson fermions in the fundamental representation of the 
gauge group. The numbers of fermion flavours were chosen according to the expected 
boundaries of the conformal window: Nf = 10 is the largest value where the theory is still 
asymptotically free, and hence it has a fixed point at small coupling where perturbation 
theory is expected to be accurate. On the other hand, Nf = 6 is expected to be close to 
the lower boundary of the conformal window, whereas Nf = 4 should be safely below it 
and hence QCD-like. 

We have measured the /3-functions of the above theories using the Schrodinger func- 
tional scheme. We also measure the mass anomalous dimension 7 for the Nf = 6 theory. 
In our analysis we have used perturbatively improved Wilson-clover action, together with 
the improvement of the boundary terms in the Schrodinger functional approach. With this 
improvement we expect most of the 0(a) errors to be eliminated. 

Overall our findings agree with the expectations. For Nf = 10 we observe a very 
small /3-function at small couplings, and a positive /3-function at larger (g 2 > 3) couplings. 
Together with the fact that due to the asymptotic freedom the /3-function must be negative 
at small enough coupling, this indicates that the /3-function has a zero somewhere at g 2 <^ 3. 
This agrees with the perturbative 2-loop /3-function which vanishes at g 2 ps 2.90. 

At the interesting case of six flavours we observe step scaling (or /3-function) behaviour 
compatible with the perturbation theory up to g 2 « 5. Above this the /3-function starts 
to approach zero, with a possible fixed point around g 2 j> 12. However, our statistical 
errors are too large in order to fully resolve the behaviour of the /3-function: instead of 
the fixed point, the /3-function could as well start to decrease again at stronger couplings. 
This kind of behaviour is characteristic for the walking coupling. It would be very in- 
teresting to resolve the behaviour in this discussed in the introduction, in the 
ladder approximation the chiral symmetry breaking sets in around the critical coupling 
g 2 ~ 47r 2 /(3C2(-R)) ~ 17. The estimated provisional fixed point is close to this value (al- 
though computed with a different schema), and thus it is indeed plausible that either the 
fixed point or the chiral symmetry breaking is realized. 

Obviously, simulation at much larger volumes (closer to the continuum) would be 
required to fully pin down the behaviour of the /3-function. However, it is just barely 
possible to reach these values of the Schrodinger functional scheme coupling with the action 
we use: in order to reach strong Schrodinger functional coupling one has to use strong bare 
lattice coupling, and the couplings used in this work are already near the largest values 
which can be used in practice. Clearly, it is important to develop an action which can be 
used at stronger couplings. 

We also measured the mass anomalous dimension using six fermion flavours. The 
dimension is seen to grow slowly as the coupling increases; somewhat slower than the 
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perturbative result at strong couplings, but overall the result in this case is as expected. 

Finally, for four flavours the observed behaviour is unsurprising: the measurements 
indicate a smoothly decreasing /3-function and QCD-like behaviour. 
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